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Abstract 



We perform a large-scale Monte Carlo simulation of the three-dimensional 
Ising model on simple cubic lattices of size L 3 with L = 128 and 256. We de- 
termine the corresponding structure factor (Fourier transform of the two-point 
function) and compare it with several approximations and with experimen- 
tal results. We also compute the turbidity as a function of the momentum 
of the incoming radiation, focusing in particular on the deviations from the 
Ornstein-Zernicke expression of Puglielli and Ford. 

PACS Numbers: 05.50.+q, 05.70.Jk, 75.10.Hk, 64.60.Fr 
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I. INTRODUCTION 



Near a phase-transition critical point, some observed quantities show a universal behavior 
that is common to a large class of systems, independently of the microscopic details. A very 
important universality class is the Ising one that is characterized by short-range interactions 
and a scalar order parameter. It describes the liquid-vapor transition in fluids, the mixing 
transition in multicomponent systems, the Curie transition in (anti)ferromagnets with axial 
anisotropy. The Ising critical behavior has been extensively studied both theoretically and 
experimentally, see Refs. ]!],§. In particular, the critical exponents, the equation of state, 
and several amplitude ratios have been determined with good precision. Another important 
quantity in the theory of critical phenomena is the static structure factor, that can be 
measured experimentally by determining the intensity of the light scattered by the fluid 
relative to the intensity of the incident light ||. To probe larger wave numbers, neutrons 
are used instead of light. At the critical density of fluids near the gas-liquid critical point or 
at the critical concentration of binary fluids near the critical mixing point, one expects for 
t = (T — T c )/T c —>■ the general scaling behavior 



s±(k) = X g±m, (i) 

where \ = C ± \t\~' y , £ is the correlation length which diverges as £ = f ± \t\~ v ', k is the 
momentum-transfer vector, and ± refers to the two phases, + (resp. — ) corresponding to 
the high- (resp. low-) temperature phase. Its absolute value is given by 

k = —sm-, (2) 

where A is the wavelength of the radiation (neutrons) in the scattering medium and 9 is the 
scattering angle. The functions g±(Q), normalized so that 

gZ\Q) = l + Cf + 0{Q*) (3) 

for Q = kC, (this defines £ as the second-moment correlation length), are universal. 
Their limiting behavior is well known. For Q small, g±{Q) is approximated by the leading 
term, the so-called Ornstein-Zernicke approximation 

go m = ^ (4) 

Such an approximation well describes the data up to Q ~ 1 and is routinely used in the 
analysis of the data with fc£ small and of the turbidity for the determination of the correlation 
length 0. On the other hand, for large Q, g±(Q) shows an anomalous decay controlled by 
the exponent r\ 

9dQ) - (5) 

Therefore, the experimental determination of the structure factor for large wavenumbers 
allows a direct determination of the exponent rj |8|-|2"0| . 
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In this paper, we compute the structure factor in the high-temperature phase for small 
values of Q by means of Monte Carlo simulations on lattices L 3 , with L = 128,256. We 
are able to determine the function g+(Q) with an error of less than 1% (resp. 2%) for 
Q ~ 5 (resp. Q < 20). These numerical results together with the most recent estimates of 
the critical exponents are then used to determine interpolations that are valid for all 
values of Q and have the correct large-Q behavior. For this purpose, we use a dispersive 
approach |2"2}-|2~3|, which allows us to determine an interpolating form for g+(Q) that agrees 
with the Monte Carlo data in the small-Q region and that well approximates (within 0.5%) 
the experimental results of Ref. |TH . 



These results are then used to compute the turbidity, i.e. the attenuation of the trans- 
mitted light intensity per unit optical path length due to the scattering with the sample. 
This quantity is routinely measured in experiments, since it allows the determination of 
the correlation length. In particular, we compute the deviations from the Puglielli-Ford 
expression [[7], that is based on the Ornstein-Zernicke approximation. 

The paper is organized as follows. In Sec. [H| we review the theoretical results for the 
structure factor. In Sec. |11 A| we define the basic observables and report the behavior 



°f 9±(Q) f° r small and large values of Q. Estimates of the constants appearing in these 
expansions are reported in Sec. |11B| . In Sec. |11 C| we discuss Bray's approximation. First, 
we discuss the high-temperature phase: we update the estimates of Ref. by using the 
most recent results for the critical exponents. Then, we generalize the approximation to the 
low-temperature phase. In Sec. |T| we discuss our high-temperature Monte Carlo results 
which are compared with approximate expressions and with the experimental data of Ref. 
T8| . In Sec. [IV] we compute the turbidity, focusing on the deviations from the Puglielli-Ford 



expression |jj due to the anomalous decay of g+(Q). We find that the turbidity is larger 
than this expression by 1% (resp. 5%) for Q = 15 (resp. 350), where Q = an d (fa is 
the momentum of the incoming radiation. 



II. THEORETICAL RESULTS 
A. Definitions 

Several theoretical results are available for the structure factor. For Q small, one can 
compute the corrections to the Ornstein-Zernicke behavior, by writing 

gg\Q) = l + Ql + Y,C±Q 2n . (6) 



n=2 



For large Q, the structure factor behaves as 



q 2 - v I Q(l-a)/v Ql/u 



a behavior predicted theoretically by Fisher and Langer |25| and proved in the field- 
theoretical framework in Refs. [BE|,p7[]. 

Beside the constants c^, the constants and defined by 
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TABLE I. Estimates of c^, S M , and S^. IHT denotes the results obtained from the analysis 
of high-temperature expansions for improved models, HT,LT results obtained from the analysis 
of high- and low-temperature expansions for the Ising model, while "e-exp." and "d = 3 g-exp." 
label the field-theoretical results, (sc) and (bcc) denote the simple cubic and the body-centered 



cubic lattice respectively. Unless stated otherwise, field-theoretical results are taken from Ref. [29], 
while the IHT estimates are taken from Ref. [21]. For we should also report the Monte Carlo 



estimate of Ref. |(J, = 0.941(11). 
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eapS 



apS ) 

si = x/(ez s 



(8) 
(9) 



are of theoretical interest. Here M gap (the mass gap of the theory) and Z gap determine the 
long-distance behavior of the two-point function in x-space: 



G{x) 



l^P_ e -M gap |x|_ 



47T \x\ 



(10) 



The critical limits of S M and S% are related to the imaginary zeroes ±iQo of g± l (Q) closest 
to the origin by 



S 



M 



si 



dg-\Q) 
dQ 2 



(11) 
(12) 



Q=±iQo 



B. Numerical results 

The coefficients c+ turn out to be very small ||, c\ ~ 1CT 4 , and this explains the success 
of the Ornstein-Zernicke approximation up to Q ~ 1. The constants c+ have been calculated 
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by field-theoretic methods. They have been computed to 0(e 3 ) in the framework of the e- 
expansion |23 , and to 0(g A ) in the framework of the d=3 g-expansion [p8fl . The perturbative 
series have been resummed in Ref. [2j| obtaining the results reported in Table |. The most 
precise estimates have been obtained from the analysis of their high-temperature expansions 
in improved models |2JJ , see the results labelled by IHT in Table |I[ 
As already observed in Ref. 



28 1, the coefficients show the pattern 



< lei", I < ... < \ct\ < 1 



for n > 3. 



(13) 



Therefore, a few terms of the expansion of g+(Q) in powers of Q 2 provide a good approxi- 
mation of g+{Q) in a relatively large region around Q = 0: as we shall see, deviations are 
less than 1% up to Q ~ 3. This is in agreement with the theoretical expectation that the 
singularity of g+(Q) nearest to the origin is the three-particle cut [^,|3|]. If this is the case, 
the convergence radius r + of the Taylor expansion of g^ 1 (Q) is r + = Since, see 

Table |, ~ 1, at least asymptotically we should have 



"n+1 



9 Cn ' 



(14) 



This behavior can be checked explicitly in the large- N limit of the TV-vector model [pgf 



The coefficients c~ are also quite small, although not as much as in the high-temperature 
case. Indeed, ~ 10~ 2 , see Table |. They have been computed using field-theoretical 
methods |J2J and from the analysis of low-temperature series [M. In the low-temperature 
phase, one also observes the pattern (|13"D, although the coefficients decrease slower. This 
is related to the fact that in the low-temperature phase the nearest singularity is the two- 
particle cut, so that convergence radius r_ of the Taylor expansion of gZ 1 {Q) is r_ = 
and therefore, 



AS 



-c n w -0.27c". 



M 



(15) 



The large-order coefficients C x , C 2 , and Cf have been computed theoretically within 
the e expansion to order e 3 |23| in the high-temperature phase and to order e 2 in the low- 



temperature phase P^j . Using the e-expansion results, we obtain 



Ct ~ 0.92, 



CI 



ci 



-2.7. 



(16) 



The corresponding low-temperature parameters C n can be derived from the high- 



temperature C n by using a set of relations derived in Ref. |27 



c; 



ct 

Ci 



(l-a)/^ 



1/1/ 



(17) 
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where 



<?.=£. 

c+ 

C-' 



Uo 



Ue = jz- (18) 

Here, C and f ± are the amplitudes of the susceptibility and of the second-moment corre- 
lation length defined above, while are defined from the critical behavior of the specific 
heat, Ch ~ ■^4 ± |^| — Using the estimates of Ref. pi| , we obtain 



C{ = 1.275(10) C+ » 1.17, 
C 2 = -0.728(5) C 2 + « -1.3, 

C 3 - = -0.345(2) C+ « 0.9. (19) 

The large-momentum behavior of the structure factor has also been studied experimentally 
and the behavior (^) has been explicitly verified in the high-temperature phase. In particular, 
the exponent rj and the constant C + have been determined. Analysis of the large-fc behavior 
of the structure factor S+(k) gives: r] = 0.017(15), C + = 0.96(4) and rj « 0.030(25), C+ » 
0.95(4) (two different parametrizations of the structure factor are used) [[J; 77 = 0.0300(15), 



C + « 0.92(1) 0; 77 = 0.042(6), C + w 0.915(21) @. No unbiased determination of C, 
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and C + is available. Fixing C + + C + = —0.9 (the e-expansion result of Ref. [p3]), Ref. [18 



obtains C + = 2.05(80) and C + = —2.95(80), in reasonable agreement with the e-expansion 
predictions. 



C. Bray's approximation 



In order to compare with the experimental data it is important to know the function 
g±(Q) for all values of Q. For the high-temperature g + (Q), several interpolations have been 
proposed with the correct large- and small-Q behavior |5| j3l|j22] -p4|,p| . The most successful 
one is due to Bray p3] , which incorporates the expected singularity structure of g + (Q). 
Here, we present Bray's interpolation together with its generalization to the low-temperature 
phase. 

In this approach, one assumes g^ 1 (Q) to be well defined in the complex Q 2 plane, with a 



cut on the negative real Q axis, starting at Q s 
r 2 = 4Sm ■ Then 



— r 



where, as discussed above, rl 



9S 



M ' 



9Z\Q) 



2 sin txt]/ 2 



± 



duu x -^F ± {u) 



S 



+ 



Q 2 



u 2 -S M u 2 + Q 2 



(20) 



where F±(u) is the spectral function, which must satisfy F±(+oo) = 1, F±(u) = for u < r±, 
and F±(u) > for u > r±. Notice the appearance of the constant Cf, which is determined, 
once F±(u) is given, by requiring g± x (0) = 1. 
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In order to obtain an approximation one must specify F±(u). Bray |23| proposed to use 



a spectral function that gives exactly the Fisher-Langer asymptotic behavior, i.e. 

P±( u ) — Q±( u ) c °t §7T77 

t~i / \n . / \ /so 



F ±.*W = p n2 , n M 2 ' 21 

P±(u) 2 + Q±{uy 



where 



P± («) = 1 + ^cos^f + ^cos— , 

^ / s Cf . Tip Ct . 7T , . 

Q±(u) = — sin — + 7-sin- , (22) 

with p = (1 — ol) J p. These definitions do not specify the spectral functions completely since 
several quantities are still unknown. First of all, we should fix the critical exponents. We 



will use the estimates of Ref. [[ZT[], obtained from the analysis of high-temperature expansions 
for improved models: 

7 = 1.2373(2), 
v = 0.63012(16), 
rj = 0.03639(15), 

a = 0.1096(5). (23) 

Several other determinations are reported in Refs. P,p9"]. For Sj^ we use the estimate 
labelled by IHT reported in Table |, while for S]^ we employ the low-temperature prediction 
of Ref. pP[| , see Table Q. We must also fix and Cjf . In the high-temperature phase, Bray 



proposes to fix C£ + Cg to its e-expansion value C£ + C% = —0.9 and then to determine 
these constants by requiring F +i #(r + ) = 0. These conditions completely fix the spectral 
function and thus the structure factor. As a check, we can compare the estimates of c+ and 
Cf obtained by using Bray's approximation g + B (Q) with the previously quoted results. We 
obtain 

C+ « 0.918, C+ « 2.56, C 3 + « -3.46, 

4 » -4.2 • 10" 4 , 4 » 1.0 ■ 10~ 5 . (24) 

The constants , C<f, and C3" are in reasonable agreement with the e-expansion results 
(PI), while and are close to the estimates reported in Table |. 

In the low-temperature phase, we have tried to follow again Bray's strategy. We have 
first set +C 3 = —0.4 and required F_ B {rJ) = 0. However, the resulting estimates of C~ 
and c~ are not in agreement with the previous results: we find Cf ~ 0.87, ~ — 1 X 10~ 3 . 
Little changes if we fix C£ + C 3 = —0.4 and use the relations For this reason, we have 
given up requiring F_ )B {rJ) = and we have simply set Ci^ = —1.3, C 3 = 0.9, as obtained 
in the previous section. Then, Bray's approximation gives 

C" « 1.0, c7 « -1.1 x 10" 2 , c7 « 1.7 x 10~ 3 , (25) 

which are close to previous estimates. A plot of Bray's approximations is given in Fig. [l|. 
Note that the structure factors in the high- and low-temperature phases are very similar. 
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FIG. 1. Scaling functions g±{Q) versus Q in Bray's approximation. We report the high- (HT) 
and low- (LT) temperature scaling functions. 

III. MONTE CARLO RESULTS 



We determine the structure factor in the region of small k — as we shall see, we are able 
to reach k ~5-10/£ by means of a large-scale Monte Carlo simulation. We consider the Ising 
model on a cubic lattice, i.e. the Hamiltonian 



H 



(Tier 



31 



(26) 



<i,j> 



where <7j = ±1 and the summation is over nearest-neighbor pairs < i, j >. We measure the 
structure factor 



S(q; P,L) = \Y, + ^ + ^ <*<PA0)<W)> 



(27) 



x,y,z 



for three different values of (3 and L: (a) L = 128, (3 = 0.2204; (b) L = 128, (3 = 0.2210; 
(c) L = 256, (3 = 0.22145. Of course, in Eq. (|27|) q = 2im/L, where n is an integer. In 



the simulation we used the Swendsen-Wang algorithm, starting from random configurations 
and discarding (2-4) x 10 4 iterations. The results of the simulations are reported in Table O. 
We report the number of iterations iV it , the susceptibility Xi the second-moment correlation 
length £ and h(q; (3, L), 



h(q-J,L)=ln 



[l + q 2 e)S{ q] (3,L) 
X 



(28) 



which directly measures the deviations from a purely Ornstein-Zernicke behavior. 

In Fig. H we plot S(q; f3,L)/x for the three lattices considered — errors are smaller than 
the size of the points — together with the experimental results of Ref. |18[] for CO2 and Bray's 
approximation. We observe good agreement, the numerical data for lattice (c) being close 
to the experimental ones. 
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TABLE II. For the three lattices considered, (a) 


, (b), and (c), we report the number 


Ul lLc! ctllUIlo 


-< v lt 


, the susceptibility Xi the second-moment correlation length £ and h(q; (3, L) for n 


- aTJOiA 

— y-w y Z7V )- 




(a) 


(b) 


(c) 


N it 


4.35 x 10 6 


3.2 x 10 6 


2.14 x 10 6 


X 


669.9(4) 


1501(2) 


6339(10) 




13.050(7) 


19.739(14) 


41.16(5) 


n 




h(q;P,L) 




1 


-0.0009(9) 


-0.0015(11) 


-0.0002(17) 


2 


-0.0002(11) 


0.0003(14) 


0.0001(25) 


3 


0.0017(12) 


0.0027(16) 


0.0019(27) 


4 


0.0039(13) 


0.0065(17) 


0.0042(27) 


5 


0.0063(13) 


0.0096(18) 


0.0067(28) 


6 


0.0093(13) 


0.0135(18) 


0.0095(28) 


7 


0.0128(13) 


0.0179(18) 


0.0123(28) 


8 


0.0178(13) 


0.0232(19) 


0.0141(28) 


9 


0.0222(14) 


0.0281(18) 


0.0179(28) 


10 


0.0270(13) 


0.0335(19) 


0.0204(28) 


11 


0.0326(14) 


0.0398(18) 


0.0234(29) 


12 


0.0383(13) 


0.0459(17) 


0.0263(28) 


13 


0.0438(13) 


0.0521(17) 


0.0290(29) 


14 


0.0510(13) 


0.0593(18) 


0.0324(29) 


15 


0.0579(13) 


0.0666(18) 


0.0353(28) 


16 


0.0647(14) 


0.0736(18) 


0.0380(28) 


17 


0.0722(13) 


0.0815(18) 


0.0409(29) 


18 


0.0806(13) 


0.0896(18) 


0.0437(28) 


19 


0.0887(14) 


0.0986(17) 


0.0478(28) 


20 


0.0975(13) 


0.1078(18) 


0.0506(29) 


21 


0.1072(14) 


0.1168(18) 


0.0538(29) 


22 


0.1158(14) 


0.1271(18) 


0.0576(28) 


23 


0.1258(14) 


0.1366(18) 


0.0616(28) 


24 


0.1367(14) 


0.1473(18) 


0.0642(29) 


25 


0.1472(14) 


0.1583(18) 


0.0676(28) 



However, at a closer look one observes tiny deviations of order 1-2%. In order to observe 
better the differences among the different approximations and data, it is useful to plot the 
function h(q; (3, L) which converges to ln[(l + Q 2 )g + (Q)] in the scaling limit. We have been 
able to observe accurately (i.e. at the level of one error bar, approximately 0.3% on g + (Q)) 
this convergence only up to Q ~ 4, as it can be seen in Fig. |3|. Indeed, only in this region we 
observe a good overlap of the results for the two lattices (b) and (c), which have the largest 
values of £. As a further check, we can compare the numerical results with the small-Q 
expansion (Q) which is expected to converge rapidly up to Q ~ 3. Using Eq. (|) to order Q 6 
(resp. Q 8 ) we obtain the curve labelled "series3" (resp. "series4") in Fig. [3[ The data (c), 
that correspond to L = 256, are in perfect agreement, confirming that in this region we are 
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FIG. 2. Function S(q;(3,L)/x versus Q = q£ for the three cases (a), (b), (c). We also report 
the experimental results of Ref. |l8[, "expt," and Bray's approximation, "Bray." 
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FIG. 3. Function h(q;f3,L) versus Q = qC for the three cases (a), (b), (c). We also report the 
experimental results of Ref. E8|, "expt," and the small-Q approximations, "series3" and "series4." 



seeing the correct asymptotic behavior. In Fig. |3] we also report the experimental results of 
Ref. | lfiH - They are systematically higher than the Monte Carlo results and indicate that, at 
least in this region, the experimental error on the structure factor is approximately of order 
0.5-1%. 

For larger values of Q, we are not able to observe scaling, as it can be seen in Fig. 
According to standard renormalization-group theory 



h(q; P, L) = h(Q, L/0 + L-"h 2 (Q, L/0 + 



(29) 



where |2T| uo = 0.83(5). Thus, we could try to extrapolate in L at L/C, fixed and then take 
the limit L/C — > oo. Lattices (b) and (c) have approximately the same L/C,, L/C ~ 6 and 
thus, in principle one should be able to extrapolate in L. In practice, corrections increase 
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FIG. 4. Function h(q;/3,L) versus Q = q£ for the three cases (a), (b), (c). We also report 
the experimental results of Ref. |fL8[ , "expt," a phenomenological interpolation , "fit," and Bray's 
approximation, "Bray." 



quickly with Q, see Fig. f|, and no reliable extrapolation can be done. In any case, we believe 
we can still use the numerical data presented in Fig. ^ to conclude conservatively that, for 
Q < 15-20, h(q; (3, L) for lattice (c) is a good approximation to the limiting function with 
an error at most of 0.02, i.e. that we can use our data (c) to compute g+{Q) with a 2% 
precision up to Q < 15-20. 

In Fig. |] we also report Bray's approximation. Such an approximation agrees nicely with 
the Monte Carlo results (c) up to Q ~ 10 and, as expected, it is lower in the region Q > 10 
where we expect the results (c) to be higher than the scaling limiting curve. Bray's function 
looks therefore a reasonable approximation to the universal scaling function, although it is 
somewhat lower than the experimental data by 1-2%. 

For the computations of the next Section, it is important to have an estimate of the 
structure factor with a reasonable error bar. For this purpose, we have determined a second 
interpolation that is in better agreement with the experimental data. We will obtain an error 
by comparing the results obtained using this interpolation and Bray's approximation. This 
interpolation may be obtained by considering expressions that agree with the numerical data 
for lattice (c) in the region Q < Q max ~ 15. We shall use again the spectral representation 
(P0|), since such an expression gives automatically the behavior ( |T4"D and ensures the correct 
small-Q behavior. In order to obtain the correct large-Q behavior, we use a generalization 
of the spectral function proposed by Bray, i.e. 

1 + J2 a - u ' n j • (3°) 

Such an expression is purely phenomenological. The first term has been introduced to 
guarantee that F &t (l) = as generally expected, while corrections of order 1/u have been 
avoided, since they would give rise to terms of order for Q — » oo that are stronger 

than those appearing in the Fisher-Langer behavior (^). In Eqs. (|20"D and (p^) we use 
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FIG. 5. Ratio t/tpf versus Qo using Bray's approximation, "Bray," and the phenomenological 
approximation, "fit." We also report the corresponding asymptotic expression r as /rpF, ("asl" and 
"as2") where r as is defined in Eq. (|34|), and the phenomenological approximation (|36|), "phen," 
valid for Q < 100. In "asl" we use C+ = 0.91797, K = 0.128735, in "as2" we use C+ = 0.92, 
K = 0.160734. 



the e-expansion estimates (O) and the values of the exponents reported in Eq. fl2"3"|). The 
constants a n are fixed by requiring 5'I 1 (0) = 1 and g+{Q) to fit the numerical data (c) up 
to Q < 15. A good fit is obtained by taking n max = 6 and a 2 = —574.128, a 3 = 7588.59, 
a 4 = —29558.9, a 5 = 43740.7, a 6 = —21715.6. The corresponding curve labelled "fit" is 
reported in Fig. [|. The results depend on Qmax used in the fit, and tend to give a lower 
curve if smaller values of Q mSuX are used. However, it is interesting to remark that, with the 
choice Q m ax = 15, the interpolation is in excellent agreement with the experimental data for 
all Q > 15, see Fig. |. 

Finally, it is interesting to remark that the Ornstein-Zernicke approximation differs at 
most 1% from the correct expression for Q < 5, while for Q > 5 the Fisher-Langer formula 



can be applied, as already observed in many experimental works, see, e.g., Refs. [ 151 , 117 ,19,20 . 



IV. TURBIDITY 



The turbidity r is defined as the attenuation of the transmitted light intensity per unit 
optical path length due to the scattering with the sample. Explicitly, it is given by 



dttS(k) 



1 - - sin 2 6 
2 



(31) 



where k = 2ko sin(6 l /2), ko = 1-nnj\ is the momentum of the incoming radiation in the 
medium, A the corresponding wavelength in vacuum, n the refractive index, and Q = (0, 9) 
the solid angle. By using Eq. (0), in the high-temperature phase we can write the turbidity 
in the form 



12 



T 



2r t- 
Ql 



2Qo 



QdQ g+(Q) 



Q 2 Q 4 



2Q 8Q 



(32) 



where Q = k C, and r is a constant that can be assumed temperature-independent in a 
neighborhood of the critical point. 

For small values of Qo, the Ornstein-Zernicke approximation can be used obtaining the 
Puglielli-Ford expression M 



tpf = T r 



2a 2 + 2a + 1 



log(2a + 1) - 



2(a + 1) 



(33) 



where a = 2Q^. 

We can also compute the behavior for large Qo by using Eq. (|7|). We obtain 

t7 2 + 2t? + 8 



2mt" 7 



C 1 +(2Q C 



77(77 + 2)(?7 + 4) 



— + K + O (q 

7] \ 



ri-(l-a)/v 




where 



/l /"OO 
QdQg+(Q) + J QdQ [ g+ {Q) - C+Q^ 2 ] 



(34) 



(35) 



In order to obtain r for all values of Qo we must use a specific form for g + (Q). We will use 
here Bray's approximation and the interpolation formula obtained using ( |3"C| ) with n max = 6, 
Qmax = 15. The difference between the results obtained using these two expressions provides 
the error on our results. In Fig. || we report r/rpp using the two different approximations 
together with their asymptotic expression r as /rpF. In Bray's approximation K = 0.128735 
while in the second one K = 0.160734. The deviations from the Puglielli-Ford behavior 



are very small and for Q > 100 are well described by the asymptotic expression (|54]) with 
C+ w 0.92 and K = 0.145(16). Estimates of the turbidity for 1 < Q < 100 can be found 
in Table |TTI|. For Q < 100 one can use the phenomenological formula 



T 



T PF 



0.666421 + 0.242399 (l + 0.0087936Q^) 018195 + 0.0911801 (l + 0.09Q^) 



. 4 \ 0.0090975 



(36) 



which is also reported in Fig. [5] ("phen"). 

We wish finally to compare our results with the approximate expressions given by Ferrell 
|34j, which require Q 3> 1 and r/logQo *C 1, i.e. 1 « Q « e 1//,? w 9 x 10 11 . By expanding 
Eq. (Q) and setting as in Ref. §L = log(4Qg) we obtain 



r 



r t~ 



Ql 



Ct{L - 1) + CU 



4 



L 

2 



K 



(37) 



In order to compare with Ferrell's results, we must compute r/(4r t 1 g{2Qo)). Since, using 
the same approximations g(2Q ) = C 1 "(2Qo) _2 (1 + i]L/2 + 0(r] 2 )), we obtain 



T 



Arot-rg(2Q ) 



r)L 2 /3 K 

- 1 — + v 7 + -pr. 1 
4 V 4 



(38) 
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rpp. We use here: (a) Bray's approximation; (b) 


general phenomenological 


interpolation based on (^ 


L]j witn n max — o ana t^ max — lo. 




Qo 


(a) 


(b) 


5 


1.004 


1.004 


10 


1.008 


1.009 


15 


1.011 


1.014 


20 


1.013 


1.017 


25 


1.015 


1.020 


30 


1.017 


1.022 


35 


1.019 


1.024 


40 


1.020 


1.026 


45 


1.022 


1.028 


50 


1.023 


1.029 


60 


1.025 


1.031 


70 


1.027 


1.034 


80 


1.029 


1.036 


90 


1.030 


1.037 


100 


1.032 


1.039 



This formula agrees with Ferrell's expression once we recognize that K — 0{rj) since K — 
for a purely Ornstein-Zernicke behavior. Numerically, we predict 3/4 + K/(r]Ci) ~ 5.1(5), 
which is smaller than Ferrell's numerical result 8.4. Ferrell's expression predicts a turbidity 
that is somewhat higher than ours. Indeed, his numerical result implies K « 0.26 in Eq. 
([34]), and as consequence we would obtain r/r PF ~ 1.06 (resp. 1.085) for Q = 100 (resp. 
1000), to be compared with our prediction r/r PF ~ 1.036(4) (resp. 1.069(3)). 

Another expression for the turbidity that takes into account the anomalous decay of the 
structure factor is given in Ref. It assumes that g+(Q) = (1 + cQ 2 )~ 1+r, l 2 , where 
c = 1/(1 - rj/2). It follows that 

[{2b + l)"/ 2 - 1] [4 - 26(77 - 4) + b\rf + 2r/ + 8)] - 4r/6(l + b) 
T = ^ ^(2 + ^(4 + ,) ' (39) 

where b = (2 — r/). Such an expression however predicts a turbidity that is too large. 
For instance, for Qo = 10 it gives r/rpp ~ 1.05, to be compared with our prediction 
r/rp F w 1.008, cf. Table |TTT]. 

Note the correct turbidity r is larger than rpp since g+(Q) decreases slower for Q — ► oo 
than the Ornstein-Zernicke approximation. However, this is apparently in contrast with the 
experimental results for the binary fluid mixture methanol-cyclohexane presented in Ref. 
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